A CLASS OF GAUSSIAN PROCESSES WITH 
FRACTIONAL SPECTRAL MEASURES 



DANIEL ALPAY, PALLE JORGENSEN, AND DAVID LEVANONY 

Abstract. We study a family of stationary increment Gaussian 
processes, indexed by time. These processes are determined by 
certain measures a (generalized spectral measures), and our focus 
here is on the case when the measure a is a singular measure. We 
characterize the processes arising from when a is in one of the 
classes of affine self-similar measures. Our analysis makes use of 
Kondratiev- white noise spaces. With the use of a priori estimates 
and the Wick calculus, we extend and sharpen (see Theorem 1 7. 1|) 
earlier computations of Ito stochastic integration developed for the 
special case of stationary increment processes having absolutely 
continuous measures. We further obtain an associated Ito formula 
(see Theorem [87T]). 
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1. Introduction 

There are two ways of looking at stochastic processes, i.e., random 
variables indexed by a continuous parameter (for example time): (i) 
One starts with a probability space, i.e., a sample space, a set Q with a 
sigma algebra B of subsets, and a probability measure P on (f2, £>), and 
a system of random variables {X(t)} on the (Q,B,P). From this one 
may then compute quantities such as means, variances, co-variances, 
moments, etc, and then derive important spectral data . These in 
turn are used in various applications, such as in solving stochastic dif- 
ferential equations. Here we are concerned with the other direction: 
(ii) Given some a priori spectral data, how do we construct a suitable 
probability space (Q, B, P) and an associated process {X(t)} such that 
the prescribed spectral data is recovered from the constructed process? 
In other words, this is a version of an inverse spectral problem. For a 
number of reasons, it is useful in the study of the inverse problem to 
focus on the case of Gaussian processes. 

A zero mean Gaussian process {X(t)} on a probability space is said 
to be stationary increment if the mean-square expectation of the in- 
crement X(t) — X(s) is a function only of the time difference t — s. 
Then there is a measure a such that the covariance function of such a 
process is of the form 

(1.1) E[X(t)X(s)*) = K a (t, s)= [ X t(u)Xs(u)*da(u), i,sel, 

JR 

where E is expectation, and we have set 

p itu i 

(1-2) Xt(u) = . 

u 

The positive measure da is called the spectral measure, and is subject 
to the restriction 

da(u) 
u 2 + 1 



;i. 3) / ^ < oo. 



The covariance function K a (t, s) can be rewritten as 



K a {t,s) = r(t)+r(s)*-r(t-s), 

2 



where 




When a is even, r is real and takes the simpler form 



, , . f 1 — cos(tu) , , . 
(1.5) r(*)= / ^^Ax(w). 



We note that some authors call spectral measure instead the measure 
u 2 da{u) rather than the measure da{u). See [28j p. 25 (7)]. 

The literature contains a number of papers dealing with these pro- 
cesses, but our treatment here goes beyond this, offering two novelties: 
the inverse problem (see above), and an operator theory of singular 
measures. Both are motivated by the need to deal with families of 
singular measures a (see (11 .11) through ( 11.41) ). Our focus is on families 
of purely singular measures a with an intrinsic spatial selfsimilarity, 
typically with Cantor support, and with fractional scaling (and Haus- 
dorff) dimension; see Section [2] below; these are measures with affine 
selfsimilarity. Note that this notion is different from self-similarity in 
the time-variable; the latter case includes fractional Brownian motion 
(fBm), studied in e.g., pfl [2J [31 HU [26] . For the latter (ffim), it is known 
that the corresponding one-parameter family of measures a consists of 
a scale of absolutely continuous measures. 

The derivative of a stationary increment process is a (possibly gener- 
alized) stationary process, with covariance function 



where a denotes the Fourier transform, possibly in the sense of distri- 
butions, of a . For a function /, recall the Fourier transform 



We note that second order stationary processes can be studied with 
the use of the theory of Hilbert spaces and of unitary one-parameter 
groups of operators in Hilbert space. One may then invoke the Stone- 
von Neumann spectral theorem, the spectral representation theorem, 
and a detailed multiplicity theory to study these processes. See for 
instance [27] . 



a(t - s) 
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An important role in the theory is played by the space Ai(a) of func- 
tions in L 2 (M, dx) such that 



This space contains in particular the Schwartz space. 

In the paper [3], see also [2], the case where a{u) is absolutely contin- 
uous with respect to Lebesgue measure, i.e., da{u) = m{u)du (where 
the Radon- Nikodym derivative m satisfies moreover some growth con- 
ditions) was considered. The study of |3] included in particular the 
case of the Brownian motion and of the fractional Brownian motion. 
A key role in that paper was played by the (in general unbounded) 
operator T m on L 2 (R, dx) defined by 



with V denoting the inverse Fourier transform, in the sense of distri- 
butions. 

In this paper we focus on the case when the spectral measure is an 
affine iterated function-system measure (AIFSs). Among the AIFS- 
measures there is a subfamily which admits an orthonormal family of 
Fourier frequencies. These are lacunary Fourier series studied first in 
a paper by one of the authors and Steen Pedersen in 1998, see [23] . 
A lacunary Fourier series is one in which there are large gaps between 
consecutive nonzero coefficients. AIFS-measures may be visualized as 
fractals in the small, while their Fourier expansions as dual fractals in 
the large. The spectral measure of such a process {X(t)} is important 
as it enters in a rigorous formulation of an associated Ito formula for 
functions f(X(t)) of the process. 

The main results of the paper may be summarized as follows: We 
construct a densely defined operator Q a from L 2 (da) into L 2 (R, dx) 
such that 



This operator Q a is the counterpart of the operator T m defined in (11.61) 
and introduced in [3]. We denote by /•—>■/ the natural isometric 
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imbedding of L 2 (M, dx) into the white noise space; See Section H] for 
details. The stochastic process {X a {t)} defined by 

X a (t) = Q^L^t]), t E R, 
has covariance function 

E[X a {t)X a {sf] = [ X t{u)Xs{u)*d<j{u) 
Jr 

Following [2E] , the measure a in this expression will be called the spec- 
tral measure of the process. Its intuitive meaning is that of "spectral 
densities", not to be confused with "power spectral measure" tradi- 
tionally used for the much more restrictive family of stochastic process, 
the stationary processes. In the case of stationary processes, and when 
the power spectral measure is absolutely continuous with respect to 
Lebesgue measure, one speaks of power spectral density (psd). It is 
then the Fourier transform of the covariance function, a function of a 
single variable, namely, the time difference. When the process can be 
differentiated, its derivative is stationary and a is absolutely continu- 
ous with respect to Lebesgue measure, its derivative is the psd of the 
derivative process. 

We show that X a (t) admits a derivative (X a (t)') d = W a (t) in the white 
noise space, which is moreover continuous in the white noise space 
norm. Furthermore 

E[W a (t)(W a (s)y] = a(t-s). 

It is found that both the processes {^(t)} and {X a (t)} are in the 
white noise space. We define a stochastic integral with respect to X a , 
and prove results similar to those of [TJ, but for a different class of 
processes studied here. 

The outline of the paper is as follows. The paper consists of nine sec- 
tions besides the introduction. The first three small sections are of 
a review nature. In Section [2] we present some material on measures 
with affine selfsimilarity. In Section [3] we recall some properties of the 
associated L 2 (da) spaces. Hida's white noise space theory is based on 
a Hilbert space, and plays an important role in our work. Its main 
features are listed in Section HJ In Sections 4-9 we develop the new 
results of the paper. In Section [5] we construct an operator for which 
( I1.7P holds. The corresponding process X a and its derivative are con- 
structed in Section [6j The associated stochastic integral and a Ito 
formula are considered in Sections [7] and |8] respectively. To contrast 
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with the measures considered here, two examples of stationary incre- 
ments Gaussian processes with measures with unbounded support are 
presented in Section [91 In Section [10] we briefly consider the case of a 
general measure a. The last section is devoted to various concluding 
remarks. 



In understanding processes {X(t)} with stationary increments, one 
must look at the variety possibilities for measures a, representing spec- 
tral measures, in the sense outlined above. Each measure-type for a 
entails properties of the associated Ito formulas for {X(t)} as it enters 
into stochastic integration formulas. 

While earlier literature on stationary-increment processes has been fo- 
cused on the case when a was assumed to be absolutely continuous 
with respect to Lebesgue measure, or perhaps the case when it is sin- 
gular but purely atomic; in this section we will focus instead on a quite 
different family of measures: purely singular and non-atomic. They 
share the following four features: 

(i) they are given by explicit recursive formulas; 

(ii) they posses an intrinsic affine selfsimilarity; see (12. ip . 

(Hi) they admit a harmonic analysis based on a lacunary Fourier ex- 
pansion; see (12 ,4p . and finally, 

(iv) the Fourier transform of a admits an explicit infinite-product for- 
mula; see ( 12. 3p . 

Definition 2.1. A Borel probability measure a on M. is said to be an 
affine iterated function system measure (AIFS) if there is a finite family 
J 7 of (usually contractive) affine transformations on R such that 



holds i.e. 

J f(x)da(x) = car ^ j- Yl J f( T ( x )) da ( x )i 
for all bounded continuous functions on K. 

The simplest examples are Bernouilli convolutions. Then card T = 2 
and there exists some fixed p > such that 



2. Measures with affine selfsimilarity 



(2.1) 




(2.2) 



T±(x) = p(x ± 1). 
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In that case, the measure da p satisfying (I2.ip . has Fourier transform of 
the form 

oo 

(2.3) a p (t) = J]cos(p fc t). 

k=l 

Note that the function 

oo 

l[cos(p k (t-s)) 

k=l 

is positive definite on the real line, since each of the functions in that 
product 

cos(p fc (t - s)) = cos(p h t) cos(p h s) + sm(p h t) sin(p fc s) 
is positive definite, and one can obtain a from Bochner's theorem. 

Cases with p of the form 

p=— , m = 2,3,4,... 
2m 

will be of special interest here. Fix m and let a m (that is with p = j-) 
be the corresponding Bernouilli measure. For t G R, set 

e t {u) = e ltu . 

Let 

( N 

(2.4) A m = 2tt I bj(2m) j , where N eN and bj G jo, y 
For instance, 

A 2 = 2tt {0,1, 4, 5, 16, 17, 20,...} 
A 3 = 27r|o,|,9,y,18,...| and 

A 4 = 2tt {0,2, 16, 18, 128, 130,...}. 
It is known that the set A m makes 

{ e A | A G A m } 

into an orthonormal basis (ONB) in Li2(da m )] we say that (a m ,A m ) is 
a spectral pair. Let us formalize this notion: 

Definition 2.2. A Borel finite measure a on R is said to have a spec- 
trum A C R if A is a discrete set and the set {e\, A G A} is an or- 
thonormal basis in L 2 (da). Then (a, A) is called a spectral pair. 
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By the Fourier basis property mentioned above, we refer to the pres- 
ence of a Fourier orthonormal basis (ONB) in the Hilbert space h 2 (da); 
and our discussion below is restricted to the case when a is assumed 
to be a finite measure. The study of these singular measures was 
initiated by one of the authors in collaboration with co-authors, see 
[2H [361 1231 [221 [2Q1 [2H UHl [51 HH El El [7]. The Fourier expansion in 
~L 2 {da) for fractal measures a differs from standard Fourier series (for 
periodic functions) in that the fractal Fourier expansion is local, much 
like wavelet expansions; see [37] for details. While the family of these 
singular measures is extensive, we found it helpful to focus our dis- 
cussion below on one of the simplest cases, the first occurring in the 
literature, see [23]. It has Hausdorff dimension= scaling dimension= 
1/2, and its support is a Cantor-subset of the real axis. 

It is proved in [23] that the support of da m is inside the closed interval 
[—1/2, 1/2], and has Lebesgue measure 0. We note however, there are 
also spectral pairs (cx m , A m ) where the measure da m is not compactly 
supported. 



For later use, we review two results from [5j El [7J [8] and [23]- Recall 
that, for i 6 K, e t {u) = e ltu . In general one does not assume that o 
has compact support. When the support of a is compact, it has a well 
defined Fourier transform, which is an entire function and not merely 
a distribution. We have the following result, proved in [221 [221 IE]- 

Theorem 3.1. Let a be a finite positive Borel measure and let A C K. 
be a discrete set. Then, (a, A) is a spectral pair if and only if 



3. The spaces L 2 (da) 



(3.1) 




Lemma 3.2. Let t, s & R. It holds that 



(3.2) 



e* - e s || L2(d(7 ) < K\t - s 



where 
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Proof: We have 

\ e t - es|| L2 (d<7) = / l e - e I d °\ u ) 

-1/2,1/2] 

|1- e lu{t - s) \ 2 da(u) 

-1/2,1/2] 



sin 



2 



4u Y~da{u] 

[-1/2,1/2] ( u{t-s) X " 



(*- S ) 



2 



< ( / M 2 rfa( M ) ) • (t - s) 2 

-1/2,1/2] 



□ 



As already mentioned, the measures a we consider are such that an 
orthonormal basis of L 2 (c?cr) is of the form 

e iXnU ,n = 0,1,..., 

where A n G 7rNo for all n e No- 

4. A BRIEF SURVEY OF WHITE NOISE SPACE ANALYSIS 

In this section we present some technical details required in the subse- 
quent sections, taken from Hida's white noise space theory. We refer 
the reader to [12], [32], P3] for more information. The facts reviewed 
here are essential for our analysis of certain stochastic integrals (Sec- 
tion [7]), and our Ito formula (Section [S]) . While convergence questions 
for stochastic integrals traditionally involve integration in probability 
spaces of paths, in our approach, the sample space will instead be a 
space of tempered distributions S' derived from a Gelfand triple con- 
struction; but there is a second powerful tool involved, a completion 
called the Kondratiev-Wick algebra. We briefly explain the justifica- 
tion for this approach below. 

The second system of duality spaces are called Kondratiev spaces, see 
Section @] for details. Further, there is a particular Kondratiev space, 
endowed with a product, the Wick product and an algebra under this 
product. It serves as a powerful tool in building stochastic integrals be- 
cause, as we show, the stochastic integral takes place in the Kondratiev- 
Wick algebra; and we can establish convergence there; see Theorem 17. II 
Moreover (see Theorem 18. ip . the stochastic integration making up our 
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Ito formula lives again in the Kondratiev-Wick algebra. 

Let S denote the Schwartz space of real-valued C°°(R) functions such 
that 

Vp, q G N , lim x p f {q \x) = 

x— >±oo 

For s G S, let ||s|| denote its L 2 (IR) norm. The function 

\\s 1 -s 2 \\ 2 

K(s 1 -s 2 ) = e 2 

is positive definite for si,s 2 running in S. By the Bochner-Minlos 
theorem (see [32], [TH Theoreme 3, p. 311]), there exists a probability 
measure P on S' such that 

(4.1) K{s) = [ e- i{s '' s) dP(s'), 

Js> 

where (s', s) denotes the duality between S and S'. Henceforth, we set 
Q = S'. The real Hilbert space W = L 2 (fi, JP, dP), where & is the 
Borelian cx-algebra, is called the white noise space. 

For s £ S and u G f2 we set 

(4.2) s{u) = (u,s) 

From (14.11) follows that the map s H- 's is an isometry from S en- 
dowed with the L 2 (M, dx) norm into W. This isometry extends to all 
of L 2 (M, dx), and we will denote the extension by the same symbol. 

We now present an orthogonal basis of W. We set £ to be the space of 
sequences («i, a 2 , . . .), whose entries are in 

N = {0,1,2,3,...}, 

where a& 7^ for only a finite number of indices k. Furthermore, we 
denote by ho, hi, . . . the Hermite polynomials. The functions 

00 

H a = H a (u) = JJh Qfc (/i fc (w)), a G t, 

k=l 

form an orthogonal base of the white noise space (the cu-dependence 
will be omitted throughout, unless specifically required). Furthermore, 
one has 

(4.3) \\H a \\ 2 w = a\, 

where we have used the multi-index notation 

a\ = «i!a 2 ! • • • , 
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The Wick product () in W is defined by the formula 



H a 0H p = H a+P , a,/3e£, 
on the basis (H a ) ae t, and is extended by linearity to W as 

(4-4) FOG = E Z^)^ 

76^ a+/3=7 

where F = J2 aee f a H a and G = J2 a€ e9aH a . See [3 Definition 2.4.1, 
p. 39]. The Wick product F()G reduces to multiplication by a constant 
when one of the elements F or G is non random. The Wick product 
is not everywhere defined in W, and one may remedy this by viewing 
W as the middle part of a Gelfand triple. The first element in the 
triple is the Kondratiev space Si of stochastic test functions, defined 
as the intersection of the Hilbert spaces T-Lk, k = 1,2,..., of series 
/ = E QG £ f<* H <* sucn that 

(4-5) 1 11/11 |I d ^$> ! )UlW a < 00. 

The third element in the Gelfand triple is the Kondratiev space 
of stochastic distributions. It is a nuclear space, and is defined as 
the inductive limit of the increasing family of Hilbert spaces T-L' k -,k = 
1,2,... of formal series Ylia&t faH a such that 

(4-6) ||/||^ d = f -^|/ a | 2 (2N)- fcQ <oo, 

where, for (3 G £, 

(2N) ±/3 = 2 ±/3l (2 x 2) ±& (2 x 3) ±/?3 
See [IU §2.3, p. 28]. 

The Wick product is stable both in S% and in S-i. Moreover, Vage's 
inequality (see [TU Proposition 3.3.2, p. 118]) makes precise the fact 
that F()G G S-i for every choice of F and G in S-\: Let I and k be 
natural numbers such that k > I + 1. Let h G %\ and u G W k . Then, 

(4.7) ||/iO«|U<A(A;-0l|/i|MHI*, 
where 

(4.8) A{k-l) = Ij2( m ) (l ~ ky 
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5. The operator Q a 

As we saw, the construction of a process {X(t)} from a fixed spectral 
measure a depends on properties of a certain operator in dx). If 

a is assumed absolutely continuous, with Radon-Nikodym derivative 
m, then this operator T m was studied earlier and it is a convolution 
operator with the square root of m. See (11.61) and [3J. In this section 
we introduce the counterpart of the operator T m in the present set- 
ting. Recall that h Q ,hi, . . . denote the Hermite functions. We define a 
unitary map W from h 2 (da) onto L 2 (R, dx) via the formula 

W(e Xn ) = h n , n = 0,l,... 

Let M u denote the operator of multiplication by the variable u in 
L 2 (do~). The formula 

(5.1) T = WM U W* 

defines a bounded self-adjoint operator from L 2 (R, dx) into itself. Fur- 
thermore, for t/> e S we set 

Q a W) = ${T)h . 

Lemma 5.1. Let ip e S. Then, it holds that: 

(5.2) (QW)(z) = Y, ( I °(y ~ x ^{y)dy] K{x) 

Proof: Using the functional calculus we have from (15.1 j) that, 

$(T) = W$(M U )W* 

and hence 

$(T)ho = W$(M u )W*h 
= W$(M U )1 
= W$(u). 

Moreover, since ip e ~L 2 (do~) we have 

00 

$( u ) = $^(^ e A„>L 2 (d<x)e An (w) 

n=0 

= £( / fWe-^^Cu)) e Xn (u), 

n=0 ^ ' 

so that W*ip is given by (15.21) . □ 
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The operator Q a is typically an unbounded operator in the Hilbert 
space L 2 (R, dx), but it is well defined on a dense domain which consists 
of the Schwartz space S C L 2 (K, dx). These facts are elaborated upon 
in the next theorem. In the statement, note that the Frechet topology 
of S is stronger than that of the L 2 (IR, dx)-norm. 

Theorem 5.2. Let ip e «S. Then, it holds that, 



(5.3) \\Q«m 2 w *)= I ivww*) 



In particular, Q a is a continuous operator from S into L 2 (R, dx). More 
precisely, 



(5.4) 



WQMWi^dx) <Vk\ { I \i>{x)\dx) + ( f \^\x)\dx 



where 



(5.5) K=l *®L <oa . 

1 + u 



Remark 5.3. For some stochastic processes it is important to relax 
condition (15.51) to 



(5.6) K p = I da( ; u i < 00 



for some p6N. In this case the estimate in (15.41) becomes 
(5.7) 

(2 2\ 1 / 2 

U\^{x)\dx^j + U\^\x)\dx 
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Proof of Theorem 15. 2t Let if> G <S. We have: 

oo „ 

\\Q^\\l 2 {R,dx) = X)| / ^(2/- A n)^(l/)0??/| 2 
n=0 

= e i yy e m(y - A "V(^w*/f 



n=0 

oo 



n=0 



E I / ( / ^)e^2/)e- mA "da(u) 

n=0 J 

oo „ 

£}| / $(u)e- iuXn da(u)\ 2 
u)\ 2 da(u), 



where we have used Fubini's theorem for the third equality, and Par- 
seval's equality for the last equality. 

We now prove ( 15. 4ft . This will prove the continuity of Q a from S 
endowed with its Frechet topology into L 2 (M, dx). For if) e S we have 



\\Q^\\l 2 (R,dx) = / \^(u)\ 2 da{u) 
Jr 



(l + « 2 )l^)| 2 ^4 
1 + 

< irmax(l + u 2 )|^(n)| 2 



<lf / \if;*if; i \(x)dx + \ \ij/*{$)'\{x)dx 
\Jr Jr 

< K ( ( I \ip{x)\dx] + ( [ \if>'(x)\dx 



where iffi(x) = ip(—x) and if)' = Furthermore, we have used 

if)*if)$(u) = \ip{u)\ 2 . 

and 

\$*ifi\{x)dx< ( I \if)\(x)dx] ( / \if$\{x)dx 



[x)dx 
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Since the expression on the right hand side in this estimate is one of 
the Frechet semi-norms of S, the continuity assertion in Theorem 15.21 
follows. The estimate f!5 .4jl further gives an exact rate of continuity. 

□ 

From equation ( 15. 3 p we can extend the domain of definition of Q G to 
a wider set, which in particular include the functions l[o,ti- This is 
explicited in the following proposition. We remark that such a result 
may be extended to more general measures ex's, see Section [TU1 

Proposition 5.4. Let f G h 2 (da) be such that, for some sequence 
(s n )neN of Schwartz functions, 

(5.8) lim |/ - tnloo = 0. 

71— >00 

Then the sequence (Q a s n )neH is a Cauchy sequence in L 2 (IR, dx). Its 
limit is the same for all sequences which satisfy (15. 8p . and will be de- 
noted by 

Q a f d = lim Q a s n . 

71— >oo 

Proof: From (15.31) follows that for every e > there is an N G N such 
that 

n,m> N \sZ. - s^|oo < e. 
Thus for such n and m 

WQaSn - QaS m \\l 2{R4x) = / ~ sZ(u)\ 2 dcx(u) 

(5 ' 9) < <r(R) ■ \s n - s^\l 

< a(R) • e 2 . 



Therefore, lim^oo Q a s n in the norm of Li2(R, cte). Call this limit qi, 
and assume that, for another sequence (t n )neN satisfying (I5.8p . we ob- 
tain another limit, say q 2 - Note that 

(5.10) lim \s n - t n |oo < lim \s n - + lim |/ - t n \^ = 0. 

n— >CO 71— >oo 71— >oo 

Then, 

\\Qi ~ Q2\\h 2 (R,dx) < \\Ql ~ QaSn\\~L 2 (R,dx) + 

+ ||Qo-Sn — Qat n \\L 2 (M.,dx) + HQcr^n ~ <?2 1| L 2 (R,dx) 
< \\qi — QaSn\\h 2 (K,dx) + 



+ a/ Ct(R) ■ \.Sn — t n |oo + \\Q a t n ~ ?2 ||l 2 (]R,<£e) , 

which goes to as n — > oo by definition of qi and q 2 and due to (I5.10p . 
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We now verify that the function 

Xt(u) = 

u 

can be approximated in the supremum norm by Schwartz functions. 
The function \t vanishes at infinity, and hence can be approximated 
in the supremum norm by continuous functions with compact support; 
see [3"51 Theorem 3.17, p. 70]. These in turn can be approximated 
by functions in S, using approximate identities, as for instance Step 
6 in the proof of Theorem 6.1 in pQ. We sketch the argument for 
completeness. Let 

1 fx 2 

(5.11) k e {x) = — =- exp ( - — 



k e is an jV(0, e 2 ) density, and therefore 

(5.12) / k e {x)dx = 1, 

Jr 

and, for every r > 



(5.13) lim / k e (x)dx = 0. 



|at|>r 



Indeed, for \x\ > r > 0, 



1 f°° _ , 1 f° x _ , 2 e 2 , 

e 2?c(z = — == / — e ^ — da; 



e f°° x 



< — == / — e 2? dx 
rV2n Jr e 2 



€ 



r 2 

-e 2? 



rv2vr 
— >■ as e -> 0. 

Theses properties express the fact that k 6 is an approximate identity. 
Applying [TO} Theorem 1.2.19, p. 25] we see that, for every continuous 
function with compact support, 

lim \\k e * f - /Hoc = 0. 
To conclude, one proves by induction on n that the n-th derivative 

is a finite sum of terms of the form 

^ exp ( - ^ p- )p(x - u)f(u)du, 



2ne Jr ^ v 2e 2 
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where p is a polynomial. All limits, 

lim x m (k e * f) (n) (x) =0 

\x\— >oo 

are then shown, using the dominated convergence theorem, and all the 
functions 

(k e *f)(x) = -^ [ eM- { ^-^)f(u)du 

V27T6 Jr t l 

are in the Schwartz space. □ 

We now compute the adjoint operator Q*. Note that it is an operator 
from L 2 (R, dx) into S', and therefore lies outside L 2 (R, c&r). We begin 
with a notation and a preliminary computation. For <fi G L 2 (IR, dx) set 



c n (<fr) = / 4>{y)K{y)dy. 

Then (c n (0)) nS N is in £ 2 . Indeed 

||(c n (0)) n6N ||| = II^IIL^)- 
We introduce the operator from L 2 (R, dx) into L 2 (da): 

CO 

(T a <f>)(u) = J2 c Me iXnU . 

n=0 

Clearly 

|H|l 2 (]R,(2x) = ||^70||L 2 (d CT )- 

Theorem 5.5. Let ip E S and <fi E L 2 (M, ofe). T/ien, 

(5.14) (Qa^,(f>)u(R,dx) = I ${u)TTM)da{u) 

J sup a 

(5.15) = / ^(y)X(0)(y))dy, 
u>/iere 

oo 

(5.16) (X(0))(y) = J^(/i„, 0) L2 (R,do ; )o : (y - A n ). 



n=0 
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Proof: We first prove (15.141) . In view of the formula (15.21) for Q a , we 

have 



(Qail), 4>)h 2 Qi,dx) = ^ yj <r{y - K)i l {y)dy S j {^J (f)(x)h n (x)dx 



f2([([ e^-^daiu)) ^{y)dy\ ( j <f>(x)h n (x)dx 

„ n \«/R \J R / / 



n=0 

oo 



f2 ( [ ( [ e iu ^(y)dy\ e-^da(u) \ ( j <j){x)h n {x)dx 
Jr \ n=0 \Jr 



{x)h n {x)dx ) ) da(u) 



ip(u) ■ T a {<j))dcr{u), 

where we have used Fubini's theorem for the third equality, and the 
continuity of the inner product for the fourth equality. 

We now turn to the second formula. The sequence ((h n , 0)L2(R,<2z))neN o 
is in i 2 . In view of ( 13. ip . the Cauchy-Schwarz inequality implies that 
X(0) in (I5.16P converges pointwise for every real y. We note that, in 
general, X(0) ^ L 2 (K, dx). For <ft and ip as in (I5.14p we have: 

(Qo4>, 0)l 2 (r,^) = j (/ °(y ~ } 'n)^{y)dy S j h n {x) \ <p(x)dx 

= V ( / a(y - \ n )ip(y)dy J ( / h n (x)4>(x)d, 

Jr \ n=0 Jr 



4>(x)h n (x)dx dy. 



To obtain the second equality, we note the following: Write 



/ a(y - X n )ip(y)dy = 
Jr Jr 



a(y - X r , 



((y 2 + l)iP(y))dy. 



y 2 + i 

Using the Cauchy-Schwarz inequality, we see that 

d( y -\ n My)dy\ 2 <d '^"^H ( I V + mmy)\' 2 dy ) . 
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In view of (13. ip . 

V) ( / S(y - \ n )i){y)dy ) h n (x) 

n=0 V* / 



belongs to L 2 (M, (ix), and we use the continuity of the scalar product. 
Furthermore we have used the dominated convergence theorem to ob- 
tain the third equality. □ 

Equations (I5.14p and (15.151) allow to compute Q* a : 

Theorem 5.6. The domain ofQ* a is the Lebesgue space L 2 (R, dx). For 
(J) e L 2 (M, dx), Qp4> is the tempered distribution defined through 



(Q*M^)s<,s= / i>{u){T a {4>))<fo{u). 



SUp <T 



Equivalently, is the tempered distribution defined by the function 

X(0) ; that is, with some abuse of notation 

oo 
n=0 



The second representation for Q* a has an important consequence: 
Theorem 5.7. It holds that 

kerQ; = {0}. 

It follows that the operator Q* a Q a is a continuous and bounded opera- 
tor from S into S'. We now provide two formulas for this operator. 

Theorem 5.8. Let ip and <p be in S. Then, 



jr(f 4>{u)e- iXnU da{u) \ ( f ^(u)e- iX ^ u da(u)\ 

n= Q \JM. / \«/R / 



(j){u)i}}{v)(kr{u). 
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Proof: By definitions of Q a (j) and of Q* we have 



= T [ W) ( [ ?(y - A„)0(l/)^ e' A "^a( 



n=0 

oo 



( / ^(n)e- ,;A ™ u (i(T(M) J ( / 



- ^n)<P{y)dy 



But 



a(y - X n )<p(y)dy = \ \J^-^da{v) ) 

${y)e- iXnV d(r{v) 



where we have used Fubini's theorem. This concludes the proof. □ 



Remark 5.9. While the operator Q a (see Theorems 15.61 through |5\8]) 
is well defined as an unbounded linear operator in the Hilbert space 
L 2 (M, dx), the other two operators Q* a and Q* a Q a are not. The reason 
is that the range of Q* is not contained in L 2 (M, dx). In fact, 

(Q* a h n )(x) = a(x - A n ), 

where a is the infinite product expression (12.31) . It can be shown that 
x i — y is not in L 2 (]R, dx); so as an L 2 (M, cfe)-operator, Q a is not 
closable (its adjoint, computed in L 2 (R, dx), does not have dense do- 
main! Nonetheless Q* in the extended sense maps L 2 (M, dx) into S'). 
The use of the ambient space S' of tempered distributions is essential. 
We illustrate the above discussion with the following diagrams: 

S L 2 (R,dx) ^-^ S' 

Qa Q* 
\ / 

S L 2 (R,dx) ^-^ S' 

In the following diagram, dom Q* is only a small subspace of L 2 (R, dx): 

L 2 (R,dx) ^-^ S' 

restricted / unbounded 

Dom Q* L 2 {R,dx) 
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6. The processes X a and W a 

In this section we build the Gaussian process with covariance function 
K a . First recall that, thanks to Proposition 15.41 the domain of the 
operator Q a has been extended to include the functions l[o,t]- We begin 
with: 

Theorem 6.1. Let Q a be as ([5]). Then, for every i,s£M, 

r ^iut g-ius y 

(Qat[0,t], QA[Q,s])h 2 (R,dx) = / da(u). 

Proof: We take (s n )neN and (i n ) n eN two sequences of elements of S 
with Fourier transforms (s^)„ e N and (i n )„eN converging in the supre- 
mum norm respectively to Xt and x* s - Then Ql[o,t] and Ql[o, s ] are the 
limit in L 2 (R, dx) of the sequences (Qs n )„ e N and (Qt n )neN respectively. 
Hence, 

(Qaho,t]iQaho,s])l,2(M.,dx) — hm (Q a S n , Qatm)L, 2 (^dx) 

n,m— >oo 

lim f Sn{u)t n da{u) 
1 e~ ius - 1 



•,iut i „—ius 



U U 



-da(u). 



□ 



Recall that we have denoted by / the natural isometric imbedding 
of L 2 (M, dx) into the white noise space. We set 

(6.1) Z n = h nj n = 0,l,... 

The Z n are independent, identically distributed J\f(0, 1) random vari- 
ables, and represented in white noise space. 

We arrive at the following decomposition: 
X a (t) = QJMt)) 



(6.2) 



Y][ °{y- K)dy ) z n . 

n=0 Vo / 



Theorem 6.2. iif ZioZds t/iat 

(6.3) \\X a {t) -X a {s)\\ w < \t-s\, t,sER. 
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The function t (->■ X a (t) is differentiable in W (white noise space), and 
its derivative is given by 



(6.4) W a (t) = J2°(t-K)Z n 

n=0 

Remark 6.3. Formulas (I6.2p and (16 .4p are analogous to, but different 
from Karhunen-Loeve expansions; see e.g., |25j . 

Proof of Theorem 16.21 : We first prove (16.31) . We have 
E a [\X a (t) - X a (s)\ 2 ] = \\X a (t) - X a (s)\\ 2 w 

oo r t 



ft 

n=0 Js 
°° ft pt 

<J2( ^dy){ \a(y - X n )\ 2 dy) 

n J S J S 



/t oo 
(J2\S(y - \ n )\ 2 )dy 
n 

(t-sf, 



n=0 



where we have used the Cauchy-Schwarz inequality and (13.1 p for the 
second and fourth equalities, respectively. 

We remark that, in view of (13. ip . W a (t) e W for every real t. We have 



X a (t) - X a (s) w _ J2n=o Il^iv ~ K)dyZ n 



t — s t — s 

oo 

- y^cr(t- \ n )dyZ n 

n=0 

Er=o {ll ( d (v ~ X n) ~ c{t - \ n )) dy} Z n 



t - s 



Er=o {f s ( e y - e t, e\ n )L 2 (d*)dy} Z n 



t - s 
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Hence, using Parseval's equality and (13. 2p . we obtain 



X a (t)-X a (s) 



t-s 



- W a (t) 



w 



< 



(t - s y 

/ \( e y- e t,ex n )L 2 (da)\ 2 dy 



n=0 



t - S 



< 



1 



< 



t-S 

(t-s) 



IK - e t \\l 2{du) dy 
(y - tfdy 



□ 



Theorem 6.4. The derivative process W a is continuous in the || • ||w 
norm. It is furthermore stationary and of constant variance, 

(6.5) E[\W a (t)\ 2 } = l. 

Proof: 

00 

\\w a (t) - w a ( s )\\ 2 w = ]T \d{t - K) - d{ s - a„)| 2 

n=0 
00 

= 5^1^* ~ e ^ e ^)-L2(da)\ 2 
n=0 

II II 2 
_ \\ e t e s\\L, 2 (da) 

= 2(1 - (e t ,e s ) h2{da) ) 
= 2(1 -d(t-s)) 

= 2 (l-f[cos 

\ n=l V 

Continuity and stationarity follow from the above chain of inequalities, 
together with the fact that {W / -(^)} is a Gaussian process. Equation 
fL5D follows from ([31]). □ 

We now turn to another type of representation for X a : 
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Theorem 6.5. Let 

H(u, u) = J2hn(oo)e iXnU G W ® h 2 (da) 



oo 



n=0 

T/ien 



g— iut i 

(X a (t))(u) = I H(u,u) da{u) 

u 



Proof: Using Fubini's theorem we have 



u 



[ [ e i{Xn - v)u da(u) 
Jo Jr 



a(X n - v). 

This together with (16.21) leads to the required conclusion. □ 



o 



7. The Wick-Ito integral 

In this section we establish a stochastic integration formula for the 
general class of stationary increment processes considered here. An 
extension of Ito's formula is addressed in section [HJ More precisely, see 
Theorem 17.11 with the use of a priori estimates and of a Wick calculus 
(from weighted symmetric Fock spaces), we extend and sharpen ear- 
lier computations of Ito-stochastic integration developed originally only 
for the special case of stationary increment processes having absolutely 
continuous spectral measures. We further obtain in the subsequent 
section an associated Ito formula (Theorem 8.2). 

The main result of this section is the following theorem, which is the 
counterpart of [U Theorem 5.1]. The fact that the derivative process 
W a (t) is W- valued (rather than lying in the larger Kondratiev space) 
allows to get a sharper statement. The proof follows the strategy of [1] 
and hence is only outlined. 

Theorem 7.1. Let Y(t), t e [a, b] be an S-i-valued function, continu- 
ous in the strong topology of S-i. Then, there exists apeN such that 
the function t i— > Y (t)0W / o -(^) is W ' -valued, and 

/ Y(t,u)0W ir (t)dt= lim ry(i kl w)o(^w)-^y. 
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where the limit is in the %' norm, with A : a = t < t\ < ■ ■ ■ < t n = b 
a partition of the interval [a,b] and |A| = maxo<fc< n -i(^fc+i — t k )- 

Proof: We proceed in a number of steps. 

STEP 1: There exists a p £ N such that Y(t) £ H' for all t £ [a, b], 
being uniformly continuous from [a,b] into H! p . 

This is proved in [Tj STEP 2 of the Proof of Theorem 5.1]. 

STEP 2: The function 1 1— > Y(t)()W a (i) is continuous over [a,b]. 

Here and in the sequel, we set || • || p =' || • to simplify notation. 
Using Vage's inequality (14.71) . it follows that, for p > 1, 

\\Y(t)OW*(t)-Y(s)QW«(s)\\ p < 

< \\(Y(t) -Y(s))0W a (t)\\ P + \\Y(8)0(W o (t) -W (8))\\ p 
<A(p)\\Y(t)-Y(s)\\ p \\W a (s)\\ + 
+ A(p)\\Y(s)\\ p \\W a (t) -W a (s)\\ 

where A{p) is defined by f)4.8p . We conclude the proof of STEP 2 by 
observing that 

II ' \\uo < II - \\W, 

which implies the continuity of the function Y{t)()W a (t) in the norm 



In view of Step 2, the integral f Y{t)()W a {t)dt makes sense as a Rie- 
mann integral of a continuous Hilbert space valued function. 

STEP 3: Let A be a partition of the interval [a, b] . We now compute 
an estimate for 

/b n-l 
Y(t)0W a (t)dt - Y(tk)0 (X a (t k+1 ) - X a (t k )) = 
k=0 

= E / (Y(t)-Y(t k ))OWAt)dt). 

As for steps 1 and 2, we closely follow pQ. Let p be as in Step 2, and 
set e > 0. Since Y is uniformly continuous on [a, b], there exists an 
T) > such that 

\t-s\ < V ^ \\Y(t)-Y(s)\\ p <e. 
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Set 

C = max \\WJs)\\o and A = A(p-N-3). 

s€[a,b] 

Let A be a partition of [a, b] with 

|A| = max{\t k+1 - t k \} < t]. 

We then have: 

J2( f tk+ \Y(t)-Y(t k ))<)WAt)dt) < 
k=o ' p 

<E / \\(Y(t)-Y(t k ))OW a (t)\\ p dt) 
< A Y,( \\(Y(t)-ntk))\\jW a (t)\\odt) 

<CA J2 r +1 \\(Y(t)-Y(t k ))\\ p dt 
< eCA{b-a), 

which completes the proof of Step 3 and the proof of the Theorem. □ 



8. An Ito formula 

We extend the classical Ito's formula to the present setting. Our present 
wider context for these stochastic processes entails important analytical 
points: Singular measures of fractal dimension, and singular operators, 
extending beyond the Hilbert space li2(M.,dx). This in turn brings to 
light new aspects of Ito calculus which we detail below. In addition 
to the examples in Section [2] (affine IFS measures), we further offer 
two examples in sect [9] below: the periodic Brownian bridge, and the 
Orenstein-Uhlenbeck processes. 

Lemma 8.1. The function 

r(t) = ||<5 CT l[0,t]llL 2 (H) 

is absolutely continuous with respect to the Lebesgue measure. 
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Proof: By (15. 3 j) we have 

ll<5al[o,t]||L 2 (K) = / \Xt(u)\ 2 da(u) 



1 — cos(tu) . , 
2 da[u) 



u- 



Since the support of da is bounded, the dominated convergence theo- 
rem allows to show that r(t) is differentiable and that its derivative is 
given by 



r'(t) = 2 I ^da(u). 
u 



Theorem 8.2. Let f : R ^ R be a C 2 (R) function. Then 
f{X c {t)) = f(X a (to)) + / f'(X o {s))0W <r {s)ds+ 



□ 



(8.1) t 

+ \ I f"{X a {s)y{s)ds, t <teR, 

1 J t 

where the equality is in the P-almost sure sense. 

Proof: We prove for t > t = 0. The proof for any other interval 
in M is essentially the same. We divide the proof into a number of 
steps. Step 1-Step 5 are constructed so as to show that (18.11) holds, 
Vt > 0, for Schwartz functions. This enables the extension to C 2 func- 
tions with compact support, with the equality holding in the 7-L' p sense. 
This implies its validity in the P-a.s. sense (actually, holding Vtu £ O), 
hence, setting the ground for the concluding step, in which the result 
is extended to hold for all C 2 functions /. 

STEP 1: For every (u,t) E R 2 , it holds that 

e %uX a {t) g yy^ 

and 

(8.2) e iuX « {t) 0WAt) e ?4 
Indeed, since X a is real, we have 

\e iuX °®\<\, Vu,teR, 

and hence e luXa ^ E W. Since W C "H' 2 , we have in particular that 
W a {t) G W 2 for all t e M, it follows from Vage's inequality (J47D that 
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(JE2D holds. 

The aim of the following two steps is formula (18 .ip for exponential 
functions. For a G R we set: 

g{x) = exp(iax). 

The proofs are as in p], taking into account that r(t) is absolutely 
continuous with respect to Lebesgue measure and are omitted. 

STEP 2: It holds that 

(8.3) g'(X a (t)) = iag(X a (t))§W*(t) + ^a) 2 g(X a (t))r'(t). 
STEP 3: Equation (18.11) holds for exponentials. 

In the following two steps, we prove (18. ip to hold for Schwartz functions. 

STEP 4: The function (u,t) i-> e mX<T ^0W / a(t) is continuous from R 2 
into H 2 . 

We first recall that the norm in W is denoted by || • See (14. 6p . The 
particular case p = in (14.61) gives in particular 

imio = Ei^i 2 - 

The structure of H' Q has been studied in [4J Section 7]. 

Recall now that the function t h- > X a (t) is continuous, and even uni- 
formly continuous, from R into W, and hence from R into H' p for any 
pGN since 

\\u\\ p < \\u\\yv, for u G W. 

The function (u,t) h-> e mX<T W is in particular continuous from R 2 into 
H 2 . Furthermore, using Vage's inequality (14.71) we have: 

|| e »^(*i)0TK(ti) -e iU2X ° (t2) 0W a {t 2 )\\ K < 

< ||( e »i^(<i) _ e ™ 2X ^' 2 ))0H^(ti)||^ + 

+ || e itt ^^)0(W ff (* 2 )-W^(ti))lk 
< A(2)||(e^^^) - e^ x ^)\\^ ■ ||^(ti)|k ; ,+ 
+ A(2)\\e^ x °W\\ n , ■ \\{W a {t 2 ) - WrihMn,, 
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where A{2) is defined by (14. 8p . This completes the proof of STEP 4 
since 

II • H2 < || • 1 1 < || • || w 

and in particular t 1— > W a (t) is continuous in the norm of T-L' and 
(u, t) 1 — y e iuX a {t) j g con tinuous in the norm of "H' 2 . 

STEP 5: (18. ip holds for f in the Schwartz space. 

The reminder of the proof is exactly as Steps 6-9 in the corresponding 
proof of [TJ Theorem 6.1], and hence omitted. □ 



9. TWO EXAMPLES 

In this section, to contrast the measures considered above, we now con- 
sider two example where the measure a has an unbounded support. 

Example (The periodic Brownian bridge). Take 

00 

(9.1) a{u) =J^5(u-2n), 

n=0 

that is the measure with support on the even positive integers with 
mass equal to 1 at each of these points. Then, it follows from the 
formula (jl.5p on r(t) that 



n=l 

Note that, 



. , *^ 1 — cos(2nt) 



\— ^ 1 — cos(2nt) 



n=l 



In view of the preceding equality, we call the associated process {X a (t)} 
the periodic Brownian bridge over [0, 7r]. We have 



'7r sin(nt) 

= V 

where, as in (16. ip . Z n = h n . 



n=l 



n 1 
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We note that, by construction, for every t, X a {t) belongs to the white 
noise space. On the other hand, the power series 



/ — oo 

' 71=1 



converges only in the Kondratiev space. More precisely, 

Theorem 9.1. Let a be given by (19. ip . For every t, we have that 
W a (t) G H' 2 and in the topology ofH' 4 

X a (t)-X a (s) 



(9.2) 



lim 

s— H 



W a {t). 



t - s 

Proof: The fact that W a (t) belongs to T-L' 2 follows from definition (14.61) 
since 

\\Z n \\\ = (2ny 2 , 

and since the Z n are mutually orthogonal in 1-L' 2 . We now turn to (19. 2p . 
We have 

X a (t)-X a (s) x /r '' 



t - s 



f*(cos(nu) — cos(nt))du 



n=0 



t-s 



But 

f*(cos(nu) — cos(nt))du f*(nu — ns)nsm(nv)du 



t - s 



t - s 
f s (u - s)du 



for some v G (s,t) 



t-s 



n \t — s\ 



and hence the limit goes to in the %' A norm since 



71 1| 4 



(2n) 



and the Z n are mutually orthogonal in W A . 



□ 



Example (The Ornstein-Uhlenbeck process). This is the solution 
of ths stochastic differential equation 

dX(t) =B{fi- X(t))dt + adB(t), t > 0, 

where B is a Brownian motion and 9 ^ 0, \x and a are parameters. We 
have 



X(t) = fi + ^W(e 



2et\-et 



29 
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and 

E[X(t)\ = e~ et + 12(1 - 
a 



E[(X(t) - E(X(t))(X(s) - E(X(s))] = ^ e -^) (e ^At _ x) 



Theorem 9.2. The centered Ornstein- Uhlenbeck process {X(t) — E[X(t)}} 
is a stationary increment Gaussian process. Furthermore 



2 

2l a ( 1 „-26»t\ 



f?[|X(t)- J B[X(t)]fl = -(l-e 



1 — cos(tw) , . , 
2 ao-(w) 



c?ct(m) 



a 2 9u 2 du 
2tt9 6 2 + u 2 ' 



10. The operator Q a for more general spectral measures 

da 

The purpose of this section is to the contrast the difference between a 
singular measure or not. This is a crucial distinction in passing from 
the given measure a to the associated process X a (t) ; i.e., in solving 
the inverse problem. The two cases are: (i) a is assumed absolutely 
continuous with respect to Lebesgue measure; versus: (ii) a is singu- 
lar. This distinction results in a dichotomy for the induced operators 
in L 2 (R, dx), so for the two operator questions. In case (i) we have a 
Radon-Nikodym derivative m , and the induced operator is T m , a self- 
adjoint convolution operator in the Hilbert space L 2 (M, dx) with the 
Schwartz space S as dense domain; see fll.6p . In the second case (ii) it 
is not, even existence is subtle. Now rather the induced operator is our 
operator Q a from Section [5j But it is much more subtle, and below we 
address some of the technical points omitted in Section [5j (a) There 
is now more than one choice for Q a ; and (b) none of the choices will 
be closable operators, referring just to the Hilbert space L 2 (M, dx). (c) 
Nonetheless, by working within the environment of the Gelfand triples 
(of Gaussian random fields), we are still able to make precise the two 
operators Q a and the corresponding adjoint Q* . For the singular case, 
i.e., case (ii), the Gelfand triple is thus essential in justifying our con- 
struction of the process {X a (t)}, existence and related properties. 



One common point between the previous work [3] and the present work 
is the construction of an operator Q a from a subspace of L 2 (1R, dx) into 
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itself (denoted by T m in |3j; see ( jl.6p ) such that 



(10.1) K a (t, s) = (Qa(l m ), Qv(Mo,s]))u(K,dx) = E[X a (t)X a ( S )*}. 



The natural isometry from L^M, dx) into the white noise space allows 
then to proceed by doing analysis in the Gelfand triple associated with 
the Kondratiev spaces, see Section [TJ Although there are numerous 
possible other Gelfand triples, Vage's inequality, see (14 .7p . is a feature 
which seems characteristic of this triple and is very useful in the com- 
putations. In the present section we present some general results on the 
existence and properties of the operator Q a . To develop the associated 
stochastic integral and Ito's formula, one has to relate the properties of 
a and of Q. Not all the arguments go through for general a's. Details 
will be presented in forthcoming publications. 

Theorem 10.1. Let a be a positive measure subject to 



and assume that dim L 2 (da) = oo. There exists a possibly unbounded 
operator Q from L2 (da) into L^M, dx), with domain containing the 
Schwartz space, such that ( 110. ID holds: 

K a (t,s) = (Q(l[0,t]),Q0-[O,s]))-L 2 QS.,dx), 



The operator Q in the preceding theorem is not unique. In the previous 
cases, a specific choice of Q was made by recipe, which depended on 
the special structure of Li2(da). In general there seems no natural way 
to chose a specific Q. We have dropped therefore the index a in the 
notation. 

Proof of Theorem I10.lt We proceed in a number of step. 

STEP 1: Let W be a unitary map from ~L 2 (da) onto L 2 (1R, dx), and let 



(10.2) 




and 



X a (t) = Ql[o,t]- 



fi = Wh, 



where b\ denotes the function 



1 



= G L 2 (da). 



VTTu' 



Then: 
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This is clear from the unitarity of W. 

STEP 2: The operator M u of multiplication by the variable u is a priori 
an unbounded operator in h 2 (da). It is densely dehned and self-adjoint 
in L12 (da). 

This follows from (15.41) . which is true for any measure a satisfying 

(HQ2D. 

It follows from the preceding step that the operator 

T = WM U W*. 
is a self-adjoint operator in Ii2(K, dx). The operator Q 
(10.3) Qi) = VT+T^^T)/! 

can therefore be computed using the spectral theorem. We claim Q 
satisfies (110. ip and (I5.3p . This is done in the next two steps. 

STEP 3: (JED holds. 

Indeed, using the spectral theorem we have 

\W\\l {W = II(v / T+W(^II 2 l M 

= [ \VTT^^(u)-=L=\ 2 da(u) 
Jr V 1 + u 

= [ mu)\ 2 da(u). 
Jr 

STEP 4: Set 

X a (t) = Ql[o,t\- 

Then 

E[X a (t)X a (s)*] = (Ql[o,t], Ql[0,s])L(K,cte)- 

Indeed, the function lp^] belongs to the domain of Q and the claim is 
a direct consequence of the isometric imbedding of L(IR, dx) inside W. 

□ 
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Corollary 10.2. Let a be a positive measure subject to 
10.4 / y-Yr<oo, 

and let {X a (t)} be the associated process as in Theorem \10.1[ Then, 
{X a (t)} has a continuous version. (By this we mean [31] that {X a (t)} 
agrees a.e. with some time- continuous process.) 

Proof: We will prove this as an application of Kolmogorov's test for 
the existence of a continuous version, [34, p. 14]. Because {X a (t)} is 
Gaussian, there exists K independent of t, s such that 

E[\X a (t) - X CT (s)| 4 ] = K (E[\X a {t) - X a (s)\ 2 ]) 2 . 
On the other hand, 

E[\X a (t) - X a (s)\ 2 ] = 2Re r(t - s) = 2 / 1 ~ ^ ~ s)u) da{u). 

Using (I10.4p we now show that 

(10.5) Re r(t) < C\t\, [t| e [0, 1] . 

The result will then follow from Kolmogorov's continuity criterion (see 
for instance [311 Theorem 2.2.3, p. 14] for the latter). 

To prove f)10.5p we proceed in a way similar as in [1] as follows. We 
compute first 

1 l-cos(fo) , , . n f 1 2 sin(f) 2 , . 

! da(u) = 2 / t 2 J 2 ' da{u) 



u 2 J t 2 U 2 

< Kit 2 (Ki > independent of t) 

<-Ki|*| for |t|G[0,l]. 

Furthermore, using the mean-value theorem for the function u t— > 
cos(tw) we have 

1 — cos(tw) = t 2 u sin(t^), £j G [0, u\. 

Thus 

1 — cos(tu) 7 / \ 9 f°° . , _ s daiu) 
^-^do{u) = t 2 I sm^—U. 



U ./i u 

da(u) 



i " Ji 

oo 

2 



< t' 

U 

< K 2 t 2 , where we use (110.41) . 
<K 2 \t\ for |t| G [0, 1], 
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where K 2 > is independent of t. Inequality (110.5j) follows and hence 
the result. □ 



In Section [5] we defined a specific operator Q a on the Schwartz func- 
tions, and then defined Ql[o t t] by approximation. Here we have used 
the spectral theorem. Still it is possible to compute Ql[o,t] via approx- 
imating sequences. We note that, in view of ( 110. 2p . the measure 

da{u) 

Mu) = ^ 

satisfies the following property: 

(10.6) Ve > 0, 3K compact and such that /i(M \ K) < e. 

The arguments in Proposition 15.41 can be adapted as follow. We take 
as a special sequence 

Sn = l[0,t] * ki/ n 

where k\/ n is defined via (15.111) . Then 

Sn{u) = xt(u) ■ e ^ . 
Instead of ( 15. 9ft we write (for the special sequence (s n )neN at hand) 

\\Q a S n - QaS m \\l 2 (R, dx ) = / \Sn{u) ~ S^(u)\ 2 d(T (tl) 

Jr 

2 2 

X t («)| 2 (e - ^ 2 ' — e~^) 2 da{u) 

r 22 
= / \Xt{u)\ 2 {e-^ - e -%*) 2 da(u) + 
Jk 

+ \Xt{u)\ 2 {e~^ - e~^) 2 da(u). 
Jr\k 

where K is a compact to be determined. 

We first focus on the second integral in the last equality above. Set 

2 2 

sup \xt(u)\ 2 (l + u 2 ){e^ - e~^) 2 = M < oo, 

m,n 
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and recall that dji(u) = . Then 

2 2 

Xt{u)\ 2 (e~~^ — e~^) 2 da(u) = 

2 2 

|Xt(w)| 2 (l + w 2 )(e~^ - e~^) 2 d/i(w) 



< M/i(R\K). 

For a preassigned e > 0, chose now the compact K such that 



fi(R \ K) < 



6 2 



2M 
Then 

..2 „,2 



: 10.7) / |^(M)| 2 (e-^ -e-^) 2 rfa( M ) < |. 



Since is compact there exists A/" such that, for n, m larger than iV, 
the integral 



„2 .,2 



e 



(10.8) / \xt(u)\ 2 (e-^ -e-^fda{u) < -. 

Therefore, for n,m > N, 

and the sequence (Q a s n )nen is a Cauchy sequence in the norm of 
L 2 (M.,dx). This provides a constructive way to compute Ql[o,t]- To 
see that the obtained limit gives the same value as the one obtained 
from the spectral theorem, it suffices to let n go to infinity in (110. 7p 
and (pTOCHD. 



11. Concluding Remarks 

We conclude with comments comparing our approach with the litera- 
ture. 



1. As in pQ, no adaptability of the integrand with respect to an un- 
derlying filtration has been made. In this sense, one may regard the 
integral defined here in fact as a Wick-Skorohod integral. 

2. Motivated in part by questions in physics, e.g., [12] and [3Q1 [33] . 
there has been a recent increase in the use of operator theory in sto- 
chastic processes, as reflected in e.g. references [TJ [2j [3j [121 1131 EH]- m 
addition, we call attention to the papers [251 EH UHl [TTJ, [29] and the 
papers cited there. In our present approach, we have been using tools 
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from the cross roads of harmonic analysis and stochastic process, as 
are covered in P [EH EE EE] ■ 
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